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Abstract 

The one-dimensional reaction diffusion process AA — > A and A$A — > 
AAA is exactly solvable through the empty interval method if the 
diffusion rate equals the coagulation rate. Independently of the par- 
ticle production rate, the model is always in the universality class of 
diffusion-annihilation. This allows us to check analytically the univer- 
sality of finite-size scaling in a non-equilibrium critical point. 
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1 Introduction 



The physics of non-equilibrium phase transitions in low dimensions is char- 
acterized by the presence of strong fluctuation effects which modify the prop- 
erties of the steady state and/or the long-time behaviour considerably with 
respect to simple kinetic equations, see |I|-pJ for recent reviews. However, 
and in contrast to equilibrium phase transition, most of the current under- 
standing of non-equilibrium phase transitions comes from numerical studies 
of certain simple models. Integrable non-equilibrium systems are still very 
rare. 

A model which has met recently with a lot of interest is the one-dimensional 
pair contact process || with single particle diffusion (PCPD). This model de- 
scribes the interactions of a single species A of particles which can undergo 
the reactions 2A — ► and AA$ —>■ AAA. If the diffusion constant d — 0, 
there is a huge number of steady states (for L sites of the order ~ <f) L , where 
4> = (1 + V / 5)/2 is the golden mean 0) and the steady state phase transi- 
tion between the active and the inactive phases is in the directed percolation 
universality class || (for d = 0, the mean pair density serves as the order 
parameter which is positive in the active phase and vanishes at the transi- 
tion towards to inactive phase). On the other hand, for d ^ 0, there remain 
just two steady states. There is a general agreement on the location of the 
transition line between the active and the inactive phases (current studies 
apply either the density matrix renormalisation group |7j or different Monte 
Carlo schemes |8],|9|]) and that the entire inactive phase should be critical and 
in the diffusion-annihilation universality class flOf . However, so far there is 



no consensus on the exponents at the active-inactive transition. 

In order to get a fresh view on this problem, we try to find a 'nearby' 
model where some analytical information might be available. In this paper, 
we shall study the following model: particles of a single species move along a 
one-dimensional lattice. Each site can be either empty (0) or else be occupied 
by a single particle (A) . Between nearest neighbour or next-nearest neighbour 
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sites the following reactions are allowed 



diffusion: 

coagulation to the left: 
coagulation to the right: 
production: 



A0 <-> 0A 
AA — > A0 
AA -> 0A 
A0A -> AAA 



rate ti 
rate d 
rate d 
rate 2d\ 



(1) 



The equality of the diffusion rate and the coagulation rates guarantees the 
integrability of the model for A = through the empty interval method [O] . 

As will be explained in the following section, the choice of this model is 
motivated by the equivalence of simple coagulation and annihilation. Con- 
sidering the above model on a finite lattice with L sites and periodic bound- 
ary conditions, we shall show in section 3 that it is integrable through the 
finite-size empty interval method []12j] . We obtain explicitly the long-time be- 
haviour of the particle density and study the finite-size scaling of the leading 
relaxation time r. We find that for all values of A, the model remains in the 
universality class of pair annihilation and we discuss the physical reasons for 
this. While that does not provide any insight into the phase transition of the 
pair contact process, the fact that A couples to an irrelevant operator allows 
us to test analytically the generalisation |H| of Privman-Fisher universality 
of finite-size scaling amplitudes []TJ] in an exactly solvable non-equilibrium 
model. Section 4 summarizes our conclusions. 



2 Motivation of the model 

Although the present understanding of the phase transition in the PCPD is 
far from being complete, it is near at hand to conjecture that the same type 
of transition occurs in a large variety of other models. More precisely, we 
expect such transitions to occur in models without parity conservation where 
(a) solitary particles diffuse, (b) particle creation requires two particles and 
(c) particle removal requires at least two particles to meet at neighboring 
sites. Examples of such reaction-diffusion processes include 

2A -> 3A, 2A -> A (2) 
2A -> AA, 2A -> A (3) 
2A -> 3A, 3A -f (4) 
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Figure 1: Typical temporal evolution of the coagulation-production process. Left: If 
particles are created to the left and to the right of a pair the system displays a phase 
transition similar to the one observed in the PCPD. Right: If particles are created between 
two particles the model is always in the inactive phase. The figure shows a typical run for 
A = 100. 

2A -> 3A, 3A -> A (5) 

A particularly interesting candidate is the coagulation/production process 
(0) . It is well known that coagulation 2 A — > A and pair annihilation 2 A — > 
are equivalent and can be related by an exact similarity transformation (for 
reviews see 0,0). Assuming that this transformation does not entirely de- 
stroy the production process 2A — ► 3A in the renormalization group sense, it 
is therefore natural to expect that the coagulation/production model exhibits 
the same type of phase transition as the PCPD. 

It is important to note that the production process 2A — > 3A in one 
dimension can be implemented in two different ways. In the standard im- 
plementation particles are created to the left and to the right of a pair of 
particles: 

• • o / o • • — > m m m . (6) 

Together with the coagulation process this implementation of the model dis- 
plays a non-equilibrium phase transition similar to the transition in the 
PCPD (see left panel of Fig. 1). In the other implementation of the pro- 
duction process, on which we will focus in the present work, a particle is 
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Figure 2: Schematic drawing of the temporal evolution of a pair of particles. For A = 
the two particles diffuse until they coagulate, forming a skeleton (bold lines). For A > 
particles are created between pairs of particles, forming a high-density region between the 
two branches. The skeleton itself remains unchanged. 

created between two other particles: 

• o« — ► • • • . (7) 

In this case the model displays a completely different behaviour. In particu- 
lar, there is no phase transition. Rather the model is always in the inactive 
phase where the asymptotic behaviour is governed by the coagulation process. 
Even the visual appearance of clusters is clearly different in both cases, as 
demonstrated in Fig. 1. Obviously, the spatial arrangement of the production 
process is crucial on a one-dimensional chain. Similar hard-core effects can 
be observed in other one-dimensional models with three-site interactions [|T^] . 

Fig. 2 illustrates why particle production between particles differs signifi- 
cantly from ordinary particle creation to the left and to the right. The figure 
sketches the temporal evolution of a pair of particles for a given realization 
of randomness. Without offspring production, the two particles diffuse until 
they meet and coagulate to a single particle. Adding the process (0), par- 
ticles can only be created between the space-time trajectories of these two 
particles. In other words, the pure coagulation process provides a 'skeleton'. 
The process (0) generates additional patches of high activity between two 
branches while the skeleton itself remains unchanged. Moreover, there is no 
way for the particles to 'cross' the branches of the skeleton. Therefore, the 
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asymptotic behaviour for t — > oo is governed by the pure coagulation process, 
i.e., we expect the density of particles to decay as t~ l l 2 . 

As can be seen in the right panel of Fig. 1, even for high values of A most 
patches with high activity die out after short time. New patches with finite 
life time can only be generated if two diffusing particles meet. As these events 
become more and more rare as time proceeds, it seems to be plausible (and 
will be proven below) that the model does not exhibit a phase transition. 

Let us point out that the PCPD with production in the middle A$A — > 
AAA and pair annihilation AA — > 00 does not have these special properties. 
Moreover, it is not exactly solvable and exhibits a phase transition similar to 
the standard PCPD. 



3 Exact solution 

We consider the model ([]]) on a chain of L sites with periodic boundary 
conditions. The exact solution can be obtained through the generalisation 



of the empty- interval method |TT],[T6[ to finite lattices |T2|]. Since we are 
interested in working out spatially averaged values of the observables, we can 
assume translation invariance from the outset. Let I n (t) be the probability 
that n consecutive sites are empty at time t. Then the mean particle density 
is given by 

p(t) = (1 - h{t)) a' 1 , (8) 
where a is the lattice constant. The equations of motion for the I n {t) are 

h{t) = 1 

^(t) = 2d[I (t)-2I 1 (t) + I 2 (t)]-2dX[I 1 (t)-2h(t)+I 3 (t)} 

^(t) = 2d[/ n _ 1 (t)-2/ n (t)+/ n+1 (t)] ; 2<n<L-l 
h{t) = (9) 

Here the boundary condition Io(t) = 1 allows one to take care of the coagu- 
lation process in the usual way, provided that the rates for coagulation and 
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diffusion coincide [|Tl|,[T2|. Obviously, the empty lattice 

hit) = 1 



(10) 



is a trivial stationary state which decouples from all other solutions. There- 
fore, we restrict our analysis to solutions with at least one particle. Since 
the last particle can never disappear, it is impossible to have L consecutive 
empty sites for a chain with L sites, leading to the other boundary condi- 
tion Iz,(t) = 0. For A = 0, eqs. (H) reduce to the known ones for simple 



coagulation with periodic boundary conditions fT2 |. 

In order to understand the A-dependent term, consider the probability 
P(nin 2 n 3 ) of realising the configuration nin 2 n% at three neighbouring sites, 
where n = • indicates an occupied and n = o an empty site. In particular, 
we have P(ooo) = I 3 and P(«oo) = P(oo«) = I 2 — I 3 . In addition, summing 
over the states of the third site 



+ P(. o o) = P(mo) 



(11) 



which yields P(« o •) = I x — 2I 2 + ^3. The production of a particle between 
two others via the process • o • only affects I\{t). 

The equations of motion are solved through the ansatz 



which leads to the eigenvalue problem 



-2dut 



(12) 
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involving an (L + 1) x (L + 1) matrix Cl. The irreversible character of the 
stochastic process reflects itself in the fact that the matrix Cl is not symmetric, 
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while probability conservation implies that the sum of the elements in a row 
of Q vanishes. Because of these properties, the real part of the eigenvalues u 
is non-negative. 

It is easy to see that the solution 



a n (0) 



n/L 



uj 







(13) 



describes the steady state with a single diffusing particle. Therefore, the 
model has two steady states, one corresponding to the empty lattice and the 
other one being the translation-invariant superposition of all single-particle 
states with an average density p av = 1/L. For the relaxational modes with 
uj > the boundary condition I (t) = 1 implies that a (uj) = 0. Similarly, 
the other boundary condition Iiif) = implies a^uj) = 0. Therefore, we go 
over to an eigenvalue problem involving an (L — 1) x (L — 1) matrix if u ^ 



-2 - A 1 + 2A 
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|T3D is solved through the ansatz 










a n {u) 




Ae ikn + Be 


-ikn . 


cu ^ 





which leads to the dispersion relation 

uj = uj(k) = 2(1 — cosfc) 



(14) 
(15) 

(16) 



and the allowed values k are obtained by inserting the ansatz ([15]) into the 
first line of eq. ( |I4"D and taking the boundary condition cll(uj) = into 
account. This leads to a system of two equations 



ik 



2e 2ik + e 61K ) + 1 



3iA: 



B(A(< 



-ifc 



2e~ 2ik + e~ 3ik 



Ae ikL + Be 



-\kL 
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which has a non-trivial solution if A; is a solution of 



tanfcl 4Asin(2fc)sin 2 (fc/2) 

tanfcL "4Acos(2A ; )sin 2 (A ; /2)-r [L7) 



We call the solutions of fllTD k m , where m — 1, . . . , L — 1. In addition, we 
can include the stationary solution by letting k = 0. Having found these, 
we can write the final result for the empty interval probabilities I n (t) in the 
form 

In(t) = ( 1 - j) + £ C m sin (& m (n - ^)) e- M ^ 4 , (18) 
V ^/ m =o 

where u m = uj{k m ) = 2(1 — cosfc m ) and the C m are real constants which must 
be determined from the initial conditions. For example, for an initially fully 
occupied lattice, one has I n (0) = <5„ i0 . If we insert the values of k m into ( Pj) 
and use (H), we obtain the average particle density as a function of time. 



Closed-form solutions of fll7|) exist for A = and A — > oo. We find 
k m = mn/L for A = and k m = mir/(L — 2) for A — > oo, respectively. For 
general values of A, we have 

m7r 27r 3 Am 3 „ , . 

fc m = -^ ^i— + ••• 5 m = 0,1,..., L-1 (19) 

Since the asymptotic scaling of the k m ~ L^ 1 is the same for A finite and for 
A = oo, even the point A = oo cannot be interpreted as being a transition 
point towards a different phase. 

From eqs. (|T8l , |l9D we see that the exact inverse leading relaxation time r 
is given by 

t~ 1 = 2du(h) ~ 2dn 2 L- 2 (l + O (l~ 2 )) . (20) 
In other words, the finite-size scaling amplitude 

A := lim L 2 t l 1 = 2dvr 2 (21) 

is independent of the particle production rate A, confirming the hand- waving 
arguments presented in Sect. 2. We point out that the value of A is equal to 
the value of the finite-size scaling amplitude of the leading relaxation time 
in the entire inactive phase in the pair contact process ||13|1 . An analogous 
universality holds for the entire spectrum of relaxation times r" 1 = 2du(k m ). 
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Figure 3: Particle density as a function of time for various values of A. All curves 
converge to a single one, demonstrating the irrelevance of the parameter A. 



Starting with a fully occupied lattice, the leading relaxation time is pro- 
portional to the time needed to reach the steady state (|i~3|) with density 
p(oo) — 1/L, see eq. (|T8|). As an immediate consequence (see section 4), the 
asymptotic decay of the particle density has to be independent of A as well. 
To verify this prediction, we performed Monte Carlo simulations. As shown 
in Fig. ^, the production process affects the curves only on a limited time 
window. Eventually all curves converge, demonstrating the universality of 
the long-time behaviour with respect to A. 



4 Conclusions 



We have seen that in our integrable coagulation-production model (g), the 
finite-size scaling of the relaxations times r m and of the steady-state particle 
density p is independent of A. We have also shown how to understand this in 
a physical way. In addition, this result can also be understood in the context 
of a recent extension Jl3| of the Privman-Fisher scaling forms fl4[] (see |L7 



for a review) to the steady states of non-equilibrium phase transitions below 
their upper critical dimension. In particular, for a ID reaction-diffusion 
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system of finite length L, the relaxation times r m should scale as (using the 
same notation as in IO) 



r m l = C L~ Z R m (o, C 2 hL 1+ *-V^) \ h=Q (22) 

and the steady-state particle density as 

p = C 2 L- p/ ^ Y' (o, C 2 hL 1+z -P /v ±) \ h _ o (23) 

where (3, u±, z are the order parameter, correlation length and dynamical 
exponents, R m and Y' are universal scaling functions, h parametrizes an ex- 
ternal source of particles, and Co, C 2 are non-universal constants ||13j| . For the 
diffusion-coagulation model at hand (when A = 0), it is known that z = 2 and 
that the time scale can be fixed by choosing Co — d | I8fl . The A-independence 



of all the r m is therefore consistent with the expected universality of the R m . 
Furthermore, the steady-state density is given by p(oo) — 1/L, which implies 
that even the generically non-universal constant C 2 does not depend on A 
and that /3/v± = 1. A simple scaling argument then shows that, in the limit 
L — > oo, also the time-dependent density 

p(t) =J2p im) e~ t/Tm ^ Po(dt)- s (24) 

where 5 = f3/v\\ = (3/(is±z) = 1/2 and the constant p is A-independent. 
That is indeed what we observe from figure [3]. 

Recall that in the pair contact process similar arguments demonstrate 
that Cq is not renormalized through the effects of the interactions, see [fL0~|, 
[X3f| and references therein. However, the universality of the amplitude A 
in the pair contact process could only be established numerically [13|]. On 
the other hand, the universality of the relaxation times in the annihilation- 
coagulation model 2A — > 0, A is trivial because of a similarity transformation 
which reduces the model to simple diffusion-annihilation. We have therefore 
obtained the first non-trivial analytic confirmation of the universality of the 
finite-size scaling amplitudes of the correlation length. 
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